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Abstract : 

The  equations  governing  the  stretching  of  a  nonlinear 
elastic  panel  of  finite  length  are  integrated  in  the  limit 
when  the  strains  are  small  but  the  strain  rates  are  comparable 
with  the  natural  frequency  of  the  panel.  It  is  shown  that 
such  a  theory  must  be  used  when  investigating  the  distortion 
and  decay  of  a  disturbance  which  is  produced  by  impact  and 
reflection.  The  mathematical  problem  involves  functions 
whose  values  repeat  at  a  sequence  of  times  which  are  separated 
by  a  time  interval  which  depends  on  the  current  value  of  the 
function. 

In  a  future  paper  we  show  that  resonant  vibrations  can 
also  be  treated  by  this  theory. 
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1.  INTRODUCTION 

Most  of  the  results  which  have  been  obtained  to  describe 
nonlinear  effects  in  wavelike  disturbances  are  valid  in  the 
limit  when  the  disturbance  is  generated  by  a  single  progressing 
wave.  These  results  are  applicable  when  the  transmitting 
medium  can  be  regarded  as  of  semi-infinite  extent  so  that  the 
effects  of  reflections  from  the  boundaries,  which  excite 
additional  components  of  the  disturbance,  may  be  neglected. 

In  this  and  future  reports  we  investigate  a  class  of  small 
amplitude  phenomena  whose  occurrence  depends  on  the  nonlinear 
interaction,  either  in  the  body  of  the  transmitting  medium 
or  at  its  boundaries,  between  different  components  of  a 
disturbance . 

Only  those  one  dimensional  disturbances  which  are 
governed  by  systems  of  totally  hyperbolic  equations  which  are 
of  the  non-dispersive  type  are  considered,  so  that  the 
equation  of  state  of  the  transmitting  medium  is  rate  independent. 
When  only  one  component  of  a  disturbance  in  such  a  system  is 
excited  it  is  generated  by  the  passage  of  a  simple  wave. 

Moreover,  according  to  classical  linear  theory  any 
disturbance  can  be  thought  of  as  generated  by  the  passage  of 
nondistorting,  nonattenuated  simple  (progressing)  waves  which, 
in  general,  interact  at  the  boundaries  but  not  in  the  body 
of  the  material.  The  mathematical  problem  consists  of 
determining  the  signals  carried  by  these  waves  from  prescribed 
initial  and  boundary  data.  It  usually  reduces  to  solving 
linear  difference  equations. 

The  small  amplitude  disturbances  which  we  investigate 
can  roughly  be  divided  into  three  classes.  In  class  I  are 
those  disturbances  for  which,  to  a  first  approximation,  there 
is  no  interaction  or  distortion  of  the  component  waves  in 
the  body  of  the  material  but,  nevertheless,  for  which  nonlinear 
effects  are  important  ?*  the  boundaries.  That  is,  although 
nonlinearity  is  of  secondary  importance  when  determining  how 
the  component  waves  propagate  it  is  of  primary  importance 
when  determining  the  signals  they  carry.  An  example  of  such 
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a  situation  occurs  when  a  column  of  gas  in  a  Kundt  tube 
is  driven  at  a  resonant  frequency  (see  Chester  [1963]). 

Class  II  disturbances  are  those  for  which  even  though 
the  components  do  not  interact  in  the  body  of  the  medium , 
nonlinearity  (amplitude  dispersion)  is  not  only  important 
at  the  boundaries  but  also  may  distort  the  signal  carried 
by  the  component  waves  so  that  shocks  may  form  and  ultimately 
dissipate  the  disturbance.  Such  phenomena  are  usually 
associated  with  the  cummulative  effect  of  locally  small 
nonlinearity.  A  typical  example  of  such  a  disturbance  occurs 
during  the  starting  problem  when  an  elastic  panelj 
initially  at  rest,  which  is  rigidly  bonded  at  one  end  is 
deformed  by  a  periodically  applied  traction  at  the  other. 
According  to  the  theory  of  linear  elasticity,  which  neglects 
all  dissipative  mechanisms,  the  resulting  motion  at  any 
particle  is  the  sum  of  an  oscillation  at  the  natural  frequency 
of  the  panel  plus  an  oscillation  at  the  frequency  of  the 
applied  traction.  In  practice,  however,  the  standing  wave 
component  is  attenuated.  In  a  future  paper,  we  show  that 
if  all  dissipative  mechanisms  are  neglected  everywhere  but 
in  shock  layers,  the  disturbance  which  occurs  is  actually 
of  class  II  and  that  our  theory  for  such  disturbances  does 
agree  with  observation,  even  at  resonance. 

Finally,  class  III  disturbances  are  those  for  which 
one  component  is  affected  by  its  interaction  with  the  other 
components  which  do  not  themselves  interact.  An  example  of 
this  is  the  Melde  phenomena  (see  Rayleigh  [1945]).  Theresa 
transverse  vibration  is  controlled  by  its  interaction  with 
a  lower  amplitude  stretching  wave.  The  frequency  of  the 
transverse  vibration  is  a  subharmonic  of  that  of  the  applied 
tension.  A  perfectly  analagous  situation  occurs  in  magneto¬ 
hydrodynamics  where  the  transverse  magnetic  field  is  controlled 
by  the  acoustic  wave.  These  disturbances  are  of  obvious 
importance  in  problems  of  control. 

In  this  paper,  as  an  example  of  a  class  II  disturbance, 
we  describe  the  evolution  of  a  small  amplitude  stretching 
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disturbance  in  an  elastic  panel  which  has  one  plane  boundary 
rigidly  bonded  when  the  other  parallel  boundary,  whose  normal 
motion  generated  the  disturbance,  is  brought  to  rest.  The 
disturbance  is  described  in  terms  of  the  variations  in  strain 
and  particle  velocity  in  the  panel  at  some  reference  time  t=0 
when  the  boundaries  are  at  rest.  Other  motions  which  are 
mathematically  equivalent  are  that  produced  in  an  inviscid  gas 
when  the  tube  in  which  the  gas  is  in  motion  is  closed  at  both 
ends,  and  that  produced  in  a  vibrating  string  after  both 
ends  are  held  fixed.  According  to  classical  linear  theory, 
if  all  dissipative  mechanisms  are  neglected,  the  resulting 
motion  is  oscillatory  with  a  period  determined  by  the  ambient 
sound  speed  a  and  width  L  of  the  transmitting  medium.  In 
practice  the  motion  is  not  exactly  periodic  and,  because  no 
energy  is  being  fed  into  the  system  at  the  boundaries, 
dissipative  mechanisms  will  ultimately  attenuate  the  disturbance. 
Which  of  the  many  possible  dissipative  mechanisms  is  mainly 
responsible  for  the  decay  will  depend  upon  the  particular 
system  as  well  as  upon  the  details  of  the  disturbance  at  t=0. 

In  part  I  of  this  report  we  describe  how  nonlinearity  (amplitude 
dispersion)  produces  large  strain  rates  and  ultimately  shocks. 

In  part  II  we  describe  how  these  shocks  attenuate  the  disturbance 
and  produce  a  uniformly  strained  panel. 

When  the  disturbance  in  the  panel  is  produced  by  an 
impulsive  motion  of  one  boundary  the  initial  disturbance  can 
be  calculated.  It  is  generated  by  the  passage  of  a  simple 
wave  (G.I.  Taylor  [1941]).  In  section  2  we  use  this  disturbance 
to  show  that  if  M  denotes  the  ratio  of  the  second  to  first 
order  elastic  constants  and  if  e(t,X)  denotes  the  strain  then, 
the  linear  theory  is  valid  only  if  the  amplitude  of  the 
disturbance  is  small  in  the  sense  that 

1 Me ! <  < 1  (1.1) 

and  the  strain  rate  is  small  in  the  sense  that 


(1.2) 
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In  the  limit  when  (1.1)  holds  but  when  (1.2)  need  not  the 
signal  distorts  and  may  generate  shocks.  In  this  small 
amplitude  finite  rate  limit^ e  is  not,  in  general,  obtained 
as  an  explicit  function  of  (X,t),  as  in  the  linear  theory: 
both  e  and  t  are  given  as  explicit  functions  of  X  and  some 
parameter  a,  which  is  invariant  at  a  characteristic  wavelet. 
Motivated  by  the  representation  of  the  solution  in  a  simple 
wave,  in  section  3  we  derive  a  representation  of  the  most 
general  disturbance  in  an  elastic  panel  in  the  small  amplitude 
limit  without  restricting  the  level  of  the  strain  rate. 

According  to  this  representation  the  disturbance  is  produced 
by  the  passage  of  two  distorting  simple  waves  which  only 
interact  at  boundaries. 

In  section  4  we  show  that  the  finite  rate  theory  must 
be  u:>M  to  discuss  the  evolution  of  the  disturbance  in  an 
elastic  panel  when  its  boundaries  are  held  fixed.  Nonlinear 
difference  equations  are  obtained  for  the  signals  carried 
by  the  component  waves.  These  are  solved  exactly  and  a 
complete  representation  is  obtained  in  terms  of  conditions 
in  the  panel  at  t=0.  As  an  application  of  the  theory,  the 
deformation  which  results  when  a  boundary  is  moved  in  an 
arbitrary  manner  over  a  time  less  than  2aQ/L,  and  then  is 
held  at  rest,  is  calculated  in  terms  of  the  prescribed 
boundary  motion  of  the  panel.  In  figures  4  and  5  we  illustrate 
how  a  disturbance  which  at  t=0  corresponds  to  a  sinusoidal 
variation  in  strain  and  zero  particle  velocity  evolves  to 
form  shocks.  In  part  II  we  show  how  these  shocks  attenuate 
the  disturbance.  We  also  generalize  our  results  to  disturbances 
which  are  governed  by  general  systems  of  hyperbolic  equations. 

2.  LIMITATIONS  OF  LINEAR  THEORY 

We  consider  one-dimensional  longitudinal  disturbance 
in  an  elastic  material  of  finite  reference  length  L.  Let 

x  =  x(X,t)  (2.1) 

be  the  position  at  time  t,  measured  from  one  boundary  X=0 
of  the  material,  of  the  material  particle  "X"  which  in  some 
reference  state,  when  the  material  was  at  constant  density  p 
and  constant  hydrostatic  pressure  pQ ,  was  at  a  distance  X 
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e  =  *,x-l.  (2 . S) 

In  particular,  for  small  strains  we  assume  that 

T  =  E[e+Me2+0(e3)] .  C2.6) 

A  special  case  of  an  elastic  material  is  an  inviscid 
polytropic  gas  which,  for  isentropic  flows,  has  an  equation 
of  state 

P  -  P0Cfi~)lr  (2.7) 

0 

where  the  density  p  is  related  to  e  by 

p(l+e)  =  pQ.  (2.8) 

According  to  (2.7)  and  (2.8),  for  an  inviscid  gas 

T  =  P0-P  =  YP^e-^^'^+CKe^)]  ,  (2.9) 

so  that 

E  =  ypo  and  M  *  -  (2.10) 
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For  algebraic  convenience  we  choose  L  as  the  unit  of 
length  and  L/aQ  as  the  unit  of  time  where 


a 


o 


(2.11) 


is  the  sound  speed  in  the  reference  state.  We  also  take  E 
as  the  unit  of  stress  and  aQ  as  the  unit  of  speed.  With 
this  scaling,  (2.3)  reads 


(e)e,x  =  u,t 


where 


(2.12) 


a2  =  T’(e)  =  l+2Me+0(e2). 


(2.13) 


Equations  (2.12)  and  (2.13)  are  supplemented  by  the  continuity 
equation 


u,x  =  e,t  (2.14) 

which  is  obtained  by  eliminating  x(X,t)  from  (2.2)  and  (2.5). 
Riemann  (see  Stoker  [1957])  showed  that  if 

c(e)  =  jSa(s)ds,  =  e{l+^Me+0 (e2) } ,  (2.15) 

then,  in  any  isentropic  wave  motion  governed  by  (2.12)  and 
(2.14),  the  combination  of  variables 

f  =  |[u-c(e)]  (2.16) 

is  invariant  at  each  a-wavelet  which  propagates  so  that 

gfl^aCe),  (2.17) 


while  the  combination  of  variables 
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g  ■=  |[u+c(e)] 


(2.18) 


is  invariant  at  each  8-wavelet  which  propagates  so  that 


DX 

Dt  |  8 


-a(e) . 


(2.19) 


If 


f  =  F(t)  at  X  =  0 • 

while  (2.20) 

g  =  G(t)  at  X  =  1, 

and  if  an  a-wavelet  is  tagged  by  the  time  if  left  X=0,  while 
a  8-wavelet  is  tagged  by  the  time  if  left  X=l,  then  according 
to  (2.16),  (2.18)  and  (2.20) 

u  =  G(8)+F(a),  and  c  =  G(8)-F(a).  (2.21) 

Classical  linear  theory  formally  appro,  mates  a(e) 
by  unity  in  (2.17)  and  (2.19).  Then,  (2.21)  imply  that 

u  "  G(8L)+F(aL),  and  c  =  G(BL)-F(aL)  (2.22) 

there  F  and  G  are  functions  of  the  forward,  and  backward 
linear  characteristic  variables 

aL  =  t-X,  and  &L  =  t+(X-l).  (2.23) 

The  statements  (2.22)  and  (2.23)  imply  that,  according  to  the 
formal  linear  theory,  the  most  general  disturbance  can  be 
regarded  as  the  superposition  of  two  components  which  do 
not  interact  in  the  body  of  the  material:  a  non-distorting , 
non-attenuated.  progressing  wave  moving  to  the  right  (the 
a -wave) ,  and  a  non-distorting  progressing  wave  moving  to  the 
left,  (the  g-wave) .  These  components  only  interact  at  the 
boundaries  of  the  material  and  the  mathematical  problem 
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reduces  to  determining  the  signal  functions  I:  and  G  from 
prescribed  boundary  and  initial  data.  This  usually  involves 
solving  linear  difference  equations. 

Since 


a(e)  =  1+Me+O (e2) ,  (2.24) 

it  might  be  thought  that  the  linear  theory  provides  a  good 
approximation  in  the  small  strain  limit  when 

| Me | <  < 1 .  (2.25) 

Then,  to  a  good  approximation,  the  traction  T(e)  and  c(e) 
can  be  calculated  from  the  strain  according  to  the  linear 
laws 


T  =  e,  and  c  =  e.  (2.26) 

However,  as  is  well  known,  the  small  amplitude  assumption 
(2.25)  is  not,  by  itself,  sufficient  to  justify  the  use  of 
linear  theory  to  determine  the  variation  in  e(X,t):  the 
strain  rates  must  also  be  "small".  The  most  relevant 
illustration  of  this  limitation,  which  is  important  in 
what  follows,  occurs  in  the  signalling,  or  starting,  problem. 
Suppose  that  the  elastic  panel  is  rigidly  bonded  at  X=1 
while  the  other  boundary  X=0  is  at  rest  for  all  times 
except  in  the  interval-  -e£t<_0  when  it  is  moved  in  a  prescribed 
way.  If  no  shocks  form,  the  pulse  which  is  generated  at 
X=0  over  this  interval  moves  towards  the  boundary  X=1  as 
a  simple  wave,  (see  G.I.  Taylor  [l^i]).  For  at  any  X, 

G(8)e0  for  times  t<t  (X)-  the  arrival  time  of  the  wavelet 
g=l-e  which  left  X=1  at  the  arrival  of  the  front  of  the 
pulse.  Then,  according  to  (2.21) 


u  =  F(ct)  and  c  =  -F(a) , 


(2.27) 
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and  (2.17)  integrates  to  give 

t  =  ct+X/a(e)  (2.28) 


for  the  arrival  time  at  X  of  the  characteristic  wavelet  a 
which  left  X=0  at  t=a.  The  wave  described  by  (2.27)  and 
(2.28)  is  amplitude  dispersed:  constant  levels  of  u  and  c 
are  carried  by  characteristic  wavelets  each  of  which  moves 
with  an  invariant  speed  which  is  determined  by  the  value  of 
c  it  carries.  Linear  theory  predicts  that  in  the  simple 
wave 


where 


c  =  *F(aL),  =  cL  say, 
t  =  oiL+X. 


\ 

► 


(2.29) 


Conditions  (2.27),  (2.28)  and  (2.29),  together  with  the 
mean  value  theorem  of  differential  calculus ,  imply  that 


Se.  ! 

F(aL)-F(cx) 

c 

c 

(aL-a) 

c 


*F  *  c  e )  I  > 


=  D (c)X 1 F ' (0)  |  , 


for  some  0 
between  'a  and  aL, 

(2.30) 


where 


D(c)  =  l-'^  l  =  |Mj  [1+0 (c) ]  .  (2.31) 

Equations  (2.25),  (2.26),  (2.30)  and  (2.31)  imply  that  the 
linear  theory,  which  neglects  amplitude  dispersion,  yields 
a  good  approximation  to  conditions  in  a  simple  wave  only  if 
the  amplitude  of  the  signal  function  is  small  in  the  sense 
that 


MF (t)  |«1, 


(2.32) 


W  Jtwvaji 
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and  if  the  rate  of  change  of  the  signal  function. is  small 
in  the  sense  that 


| MF *  (t)  j  «1.  (2.33) 

Condition  (2.32)  states  that  the  amplitude  of  imposed  strain 
must  be  small  compared  with  the  magnitude  of  the  ratio  of 
first  to  second  order  elastic  constants,  (=  for  a  gas). 
Condition  (2.33)  states  that  the  amplitude  oi  the  product  of 
travel  time  of  the  pulse  front  and  the  imposed  strain  rate 
must  also  be  small  compared  with  the  magnitude  of  the  ratio 
of  first  to  second  order  elastic  constants.  Theories  for 
which  the  restriction  (2.32)  holds  but  the  restriction  (2.33) 
does  not  are  called  small  amplitude ,  finite  rate  theories . 

Such  theories  were  first  used  to  solve  significant,  problems 
in  gasdynamics  by  Whitham  (|<?S4>)  .  They  have  been  generalized, 
and  used  to  solve  problems  governed  by  general  systems  of 
hyperbolic  equations  by  Varley  and  his  co-tirorkers  0<?bS)  , 

(Kfi>l)  and  (l^l^)  .  Up  to  now  these  theories  have  only  been 
used  to  investigate  disturbances  which  are  generated  by 
progressing  waves,  as  in  the  signaling  problem.  In  this 

report  we  show  that  such  theories  can  be  generalized  to  account 
for  effects  due  to  reflected  waves.  As  an  application, 
we  show  that  such  a  theory  is  necessary  when  investigating  the 
decay  of  a  standing  wave,  or  the  "shut-off"  problem  in  an 
elastic  panel.  In  future  reports  we  will  show  that  resonant 
wave  motions  can  also  be  analyzed  by  the  theory. 

For  future  reference,  and  to  motivate  the  mathematical 
approach  we  use  to  generalize  the  theory,  we  briefly  review 
conditions  in  a  small  amplitude  but  finite  rate  simple  wave. 

To  a  first  approximation,  according  to  (2.26)  and  (2.27) 

c  =  e  =  -u  =  -F(a)  (2.34) 


where  the  arrival  time  t=t(a,X)  of  the  characteristic  wavelet 


-li¬ 


ce  is  given  by  (2.28),  and  (2.24),  as 

t-X  =  ct+MF(cs)X.  (2.35) 

Because  the  a -wave  described  by  (2.34)  and  (2.35)  is  amplitude 
dispersed,  its  profile  distorts  as  it  propagates.  A  measure 
of  this  distortion  is  the  incremental  arrival  time  of  the 
a-wavelet  at  X, 


p(<x,X)  =  t,a  =  aL>a  =  1+MF '  (a)X.  (2.36) 

At  any  a-wavelet,  the  ratio  of  the  strain  rate  to  input 
strain  rate 


e  ’  t^e  ’t | X-0  e,x/e,x|x=0  p  *  (2.37) 

According  to  linear  theory  p=l  and  the  level  of  the  strain 
rate  in  the  material  is  always  bounded  by  its  value  at  input. 
This  theory  gives  a  valid  approximation  when 


X«|MF'  (a)  I"1. 


(2.38) 


By  contrast  the  finite  rate  theory  predicts  that  if  | F ' (a) | > |M| -1 
at  some  wavelet  then,  before  the  a-wave  reaches  X=l,  p-*0 
at  some  (X,t)  and  the  profile  distorts  so  much  that  the  strain 
rate  becomes  unbounded  compared  with  its  level  at  input  and 
at  least  one  weak  shock  forms.  The  main  effect  of  such  shocks 
is  to  attenuate  the  disturbance.  We  will  show  in  part  II 
that  in  an  elastic  material,  where  all . attenuating  mechanisms 
are  neglected  everywhere  but  in  shock  layers,  it  is  this 
shock  attenuation  which  is  actually  responsible  for  the  decay 
of  a  standing  wave. 

Shocks  attenuate  a  simple  wave  because  the  characteristic 
wavelets,  each  of  which  carries  a  constant  value  of  e  and  u, 
coalesce  into  it.  In  fact,  according  to  the  elastic  model, 
if  a  shock  is  allowed  to  propagate  into  an  "infinite  region" 
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all  wavelets  except  those  carrying  vanishingly  small  values 
of  u  and  e  will  coalesce  into  it,  or  some  other  shock,  so 
that,  ultimately  the  disturbance  will  be  fully  attenuated.  The 
equations  governing  conditions  at  a  shock  in  an  elastic 
material  are  well  known  (see  Varley  \Vtis"0) .  Because  they  are 
needed  in  what  follows,  we  briefly  indicate  how  they  are 
derived.  If  t=W(X)  denotes  the  arrival  time  of  the  shock 
at  X,  if  <x+(X)  and  a~(X)  are  the  characteristic  wavelets 
immediately  ahead  and  behind  the  shock  at  X,  then  the  usual 
jump  conditions  at  a  weak  shock  imply  that 

W*(X)  =  1  +  ^M{F(ct+)+F(cT)}  (2.39) 

while  the  condition  that  a+  and  a”  are  at  the  shock  at.  the 
same  time,  together  with  the  statement  (2.34),  also  imply 
that 


W-X  =  a++MF (a+)X  =  cT+MF(a')X.  (2.40) 


Once  the  input  signal  F(a)  is  specified,  equations  (2.39) 
and  (2.40)  govern  the  variations  in  W(X),  a+(X)  and  a“(X) 
at  any  shock.  For  the  special  case  of  a  shock  moving  into 
an  undisturbed  region,  where  F(a  )h0  for  a  <-e,  these  equations 
integrate  to  give 


X 


2 

MF2(cT) 


■a 

F(a)da, 

-e 


and  t 


X+a~+MF(a“)X 

(2.41) 


as  the  shock  trajectory.  According  to.  (2.41),  only  those 
w avelets  at  which  MF<0  can  coalesce  into  this  shock.  If 
a=0  is  the  next  zero  of  F(a)  after  a=-e,  then  the  amplitude 
of  the  pulse  generated  at  X=0  over  the  time  interval  -e<t<0 
is  vanishingly  small  at  distances  which  are  large  compared 
with  the  attenuation  length, 
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SLa  =  -2M  F(s)ds ,  (2.42) 

*  -e 

of  the  pulse.  For,  according  to  (2.41),  as  X/2,  •*■«>  the  amplitude 
of  this  pulse,  and  the  shock, 

-MF(cf)*(*  /X)1/2.  (2.43) 

Note  that  the  decay  in  strength  of  the  pulse  as  it  propagates 
is  determined  by  SL  which  is  a  measure  of  the  total  energy 
content  of  the  input  signal. 

3-  SMALL  AMPLITUDE,  FINITE  RATE  THEORY. 

In  this  section  we  show  that  when  both  the  signal  functions 
F  and  G  are  not  identically  zero,  the  linear  approximation, 

(2.22)  and  (2.23),  is  valid  only  when  their  amplitudes  are 
small  in  the  sense  that  both 

|MF (t)  | «1 ,  and  |MG(t)|«l  (3.1) 

while  their  rates  of  change  are  small  in  the  sense  that  both 

1 MF *  (t)  |  «1 ,  and  |MG'(t)|«l.  (3.2) 

In  (3.2)  a  dash  denotes  differentiation  with  respect  to  time 
t,  whose  scale  of  measurement  has  been  chosen  so  that  the 
travel  time  of  a  sonic  front  from  boundary  to  boundary  of  the 
material  is  unity.  We  integrate  equations  (£*12.  )  and  (2-14-) 
in  the  small  amplitude  finite  rate  limit  when  the  restrictions 
(3.1)  hold  but  when  (3.2)  need  not.  In  contrast  to  the 
linear  theory,  the  finite  rate  theory  predicts  that  amplitude 
dispersion  distorts,  and  that  weak  shocks  attenuate  the 
disturbance.  However,  the  two  theories  have  one  important 
feature  in  common:  the  a-wave  and  $-wave  do  not  interact  in 
the  body  of  the  material.  It  is  this  one  single  feature  which 
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makes  the  standing  wave  problem,  as  well  as  that  of 
resonant  motions,  mathematically  tractable.  Even  though, 
in  general,  the  problem  of  determining  the  signal  functions 
F  and  G  from  the  prescribed  initial  and  boundary  data 
involves  solving  nonlinear  difference  equations,  it  can 
be  solved  for  these  problems  to  show  that  the  assumptions 
(3.2),  inherent  in  the  linear  approximation,  are  invalid 
and  that  the  finite  rate  theory  must  be  used. 

Since  the  equations  governing  any  disturbance  are 
linear,  with  known  coefficients,  if  the  hodograph  variables 
(f.g^and  (t,X)j  are  used  as  independent,  and  dependent* 
variables,  this  choice  of  variables  might  seem  the  natural 
one.  However,  the  boundary  value  problem  associated  with 
the  class  of  disturbances  considered  here,  which  requires 
that 


when  X  =  0,  t  =  F-1(f) 

while 

when  X  =  1,  t  =  G-1(g) 

is  a  formidable  free  boundary  value  problem  in  terms  of 
these  variables  even  when  the  signal  functions  F(t)  and 
G(t)  are  known.  The  hodograph  description  also  has  two 
other  disadvantages.  It  is  difficult  to  interpret  any 
physical  approximation  in  terms  of  these  variables.  For 
example,  the  disturbance  generated  by  a  simple  wave,  which 
is  characteristic  of  the  disturbances  studied  here,  cannot 
be  described  in  a  straightforward  way  in  the  hodograph 
formalism.  The  second  disadvantage  is  that  the  use  of  the 
hodograph  technique  is  limited  to  disturbances  which  are 
governed  by  a  system  of  two  first  order  hyperbolic  equations. 

The  procedure  used  here,  which  can  be  extended  to 
general  systems  of  hyperbolic  equations,  is  to  choose 
variables  which  are  optimal  for  the  particular  component 
of  the  disturbance  which  is  being  investigated.  These 
variables  are  different  for  the  a-wave  component  and 
*1* 

Ludford  [1952]  tackled  this  formidable  problem.  He  used 
the  hodograph  approach  to  calculate  the  time  at  which  shocks  form. 
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8-wave  component  of  the  disturbance.  This  approach  works 
because,  to  a  first  approximation,  these  components  do  not 
interact  in  the  body  of  the  material. 

3.1  The  a-wave. 

Motivated  by  the  representation  (2.34)  and  (2.35), 
to  describe  conditions  in  the  a-wave  we  use  (a,X)  as  basic 
independent  variables.  As  basic  dependent  variables  we 
take  the  incremental  arrival  time 

p(a,X)  =  t,a(u,X),  (3.3) 

and  c(a,X).  In  terms  of  these  variables 

u  =  c+2F (a) .  (3.4) 

Differentiating  (2.17)  with  respect  to  a  at  constant  X, 
and  considering  a(c)  as  a  function  of  c  rather  than  e, 
yields 


P’x  *  c>«  ■  °-  C3.S) 

X  a2(c)  “ 

Differentiating  the  last  of  equations  (2.21)  with  respect  to 
a  at  constant  (3  yields 

Da|g  =  C ’ a  +  Da|B  C,X  = 


Since 


Dt 

na|e 


P+a_1  n||3  from  (2.17), 

'a_1  Sale  £rom  (2,19^ » 

DX  a 

Da|e  "'IP* 


(3.7) 


(3.8) 
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If  (3,8)  and  (5.6)  are  used  to  eliminate  c,^,  condition  (3.5) 
yields 

_3_ 

(/aTc)  p),x  =  [a(c)]“2a' (c)F' (a) ,  (3.9) 

which  is  regarded  as  an  equation  governing  the  variation  of 
p  at  each  a-wavelet.  The  characteristic  condition  (3.9) 
is  supplemented  by  (3.6)  which,  using  (3.8),  reads 

§f | 3  =  c,a  -  |a(c)pc,x  =  -F  '  (a) .  (3.10) 

Equations  (3.9)  and  (3.10),  which  are  two  first  order 
equations  for  p(a,X)  and  c(a,X) ,  are  in  a  form  which  is 
specially  suited  for  investigating  conditions  in  an  a-wave 
and,  in  particular,  how  it  is  affected  by  its  interaction 
with  the  8-wave.  Suppose,  for  example,  we  wish  to  determine 
conditions  in  the  a-wave  which  is  generated  at  X=0  for  times 
t>0.  Conditions  in  this  wave  are,  in  theory,  determined  by 
the  values  of  F(a)  and  c(a,l)  at  values  of  a_>0  together  with 
the  value  of  c(0,X)  for  0<X_<  1 ,  (see  fig.  1).  We  say  that 
the  interaction  of  the  8-wave  with  the  a-wave  can  be  neglected 
if  the  trajectory  of  any  wavelet  is  determined  by  F(a) 

for  and  is  independent  of  c(a,l)  and  c(0,X).  There 

is  no  such  interaction  in  the  small  amplitude,  finite  rate 
limit.  For,  without  restricting  F’(a),  integrating  (3.9) 
from  X=0  and  using  the  mean  value  theorem  of  integral 
calculus  together  with  the  restriction  | Me | <<1  to  estimate 
the  error  yields 


p  =  1+MF 1 (a)X  (3.11) 

as  a  uniformly  valid  first  approximation  for  the  incremental 
arrival  time  p=t,  .  Equation  (3.9)  and  the  mean  value 
theorem  also  implies  that  (3.11)  can  be  formally  integrated 
to  give 
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t  =  X+a+MF  fct)  X  (3.12) 

which  is  a  uniformly  valid  implicit  equation  for  the  first 
approximation  to  a(t,X).  According  to  (3.12),  just  as  in  a 
simple  wave,  the  trajectory  of  any  a-wavelet  depends  only 
on  the  signal  F(a)  it  carries. 

To  obtain  higher  order  approximations  for  conditions 
in  the  a-wave  the  iteration  scheme 

3 

(/aTcT)  Pn),x  =  [a(cn)]-2a’ (cn)F’ (a)  (3.13) 

and 

=  t3-14) 

with 

cQ  =  0,  and  pQ  =  l+MF’(a)X  (3.15) 

may  be  used.  Higher  order  approximations  may  also  be 
obtained  as  the  limiting  case  as  |MF(ct)|-»-0  of  the  regular 
perturbation  scheme  in  which  the  Riemann  invariant 

r  =  c+F(a)  =  e[ro(a,X)+er1(a,X)+0(e2)]  (3.16) 

and 


r  Po(a,x)  +  ep1(a,X)+0(e2) 


(3.17) 


where 


e  =  max . | r | . 

Since  in  a  simple  wave  r=0,  (3.16)  and  (3.17)  describe  a 
perturbed  simple  wave.  Note  that  the  perturbation  is  only 
regular  if  r  and  p  are  regarded  as  functions  of  X  and  the 
exact  characteristic  variable  a. 
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When  an  analysis  similar  to  that  described  by  equations 
(3.j)-(3.12)  Is  applied  to  the  8-component  of  the  disturbance 
using  (8,X)  as  independent,  and  (t,c)  as  dependent,  variables ? 
a.  complete  representation  of  conditions  in  any  small  amplitude, 
finite  rate  disturbance  is  obtained.  It  is  given  by 

u  =  G(0)+F(d)  (3.18) 


and 


c  =  G(8)-F(ct)  , 

(3.19) 

where  a(t,X)  and  g(t,X)  are  given  implicitly  by 

t-X  =  d+MF ( a) X 

(3.20) 

and 

t+(X-l)  =  8+MG(a)  (X-l)  . 

(3.21) 

If  a(t,X)  is  computed  from  (3.20)  then,  in  the  limit 

when  MF(d)-*-0  but  MF'(a)  is  unrestricted,  F(a(t,X))  is  a  uniformly 

valid  first  approximation  for  F  as  a  function  of  (t,X).  More- 

3t  3 1 

over,  to  obtain  the  first  approximation  to  ^  and  in  this 
limit,  equation  (3.20)  may  be  formally  differentiated  to  give 


lim 

at 

=  lim 

1+MF  (d)  =  1 

MF-+0 

MF->-0 

lim 

ii 

3d 

=  lim 

[1+MF’ (d)X] 

MF-+0 

MF-»-0 

'  =  1+MF 1 (a)X. 

To  a  second  approximation,  of  course, 

|t  =  2. -Me  =  1+M[F (a)  -G(g)]  . 

The  linear  theory  formally  approximates  a  and  8  in  (3.20) 
and  (3.21)  by  the  explicit  expressions 

aT  =  t-X  (3.22) 

ii 

and 

8l  =  t+(X-l)  (3.23) 
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Thi's  is  permissable  i£  and  only  if  conditions  (3.1)  and  (3.2) 
are  satisfied. 

4.  DISTORTION  OF  A  STANDING  WAVE 

In  this  section  we  use  the  representation  (3 . 18) - (3 . 21) 
to  describe  the  disturbance  in  an  elastic  panel  after  its 
boundaries  X=0  and  X=0.  are  brought  to  rest.  This  disturbance  is 
described  for  all  t>0  in  terras  of  the  variations  in  <J>(X)  and 
tJ»(X),  0<X<_1,  the  values  of  F  and  G  in  the  panel  at  t=0.  These 
are  related  to  u  and  c  in  the  panel  at  t=0  by  (3.18)  and  (3.19). 
The  nonlinear  mixed  initial  and  boundary  value  problem  which 
is  solved  here  reduces  to  determining  the  signal  functions  F(a) 
and  G(8)  from  (3.20)  and  (3.21)  in  the  region  0<X<1  for  all 
t>_0  from  the  initial  data  that  when  t=0, 

F  =  4> (X)  and  G  =  i|»(X)  ,  0<X<1,  (4.1) 

and  from  the  boundary  data  that  on  X=0  and  X=1 

u  =  G (8)+F (a)  =  0,  all  t>0.  (4.2) 

For  comparison,  and  to  motivate  the  steps  we  take  to  solve 
the  nonlinear  problem,  we  review  the  argument  used  to  construct 
the  solution  to  the  linear  problem.  Reference  to  fig.  2  will 
be  helpful. 

4 . 1  Linear  theory 

Linear  theory,  according  to  (3.22)  and  (3.23),  predicts 
that  in  the  region  0<X<1  when  t>0  the  characteristic  variables 
(a, 8)  range  from  -1  to  “.  To  determine  the  signal  functions 
F(a)  and  G(8)  over  this  range,  note  that  the  condition  that 
(4.2)  holds  on  X=0  when,  by  (3.23), 

8  =  a-1,  a>T),  8^-1  (4.3) 

implies  that  F  and  G  are  related  by  the  condition  that 

F(R)+G(R-1)  =  0  when  R> 0 .  (4.4) 

Similarly,  the  condition  that  (4.2)  holds  on  X=1  when,  by 
(3.22)  , 
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o  =  8-1  (4.S) 

implies  that 

F(S-1)+G(S)  =  0  when  S>0.  (4.6) 

I£  G  is  eliminated,  by  taking 

R  =  S+l,  (4.7) 

(4.4)  and  (4.6)  imply  that 

F (S+l)  =  F(S-l)  when  S  >  0 ,  (4.8) 

or,  equivalently,  that 

F  (a+2)  =  F (a)  when  a>-l.  (4.9) 

Condition  (4.9)  states  that  F(a)  is  periodic  with  period 
2  for  all  cp^-1:  if 

F(a)  =  f(a)  when  -L<a<l  (4.10) 

then 

F(a)  =  f(a-2m)  when  2m-l<a<2m+l ,  m=0,l,2,..  .  (4.11) 

Similarly,  G (3)  is  periodic  with  period  2  for  all  8_>-l  so 
that  if 

G(8)  =  g(8)  when  -1<8<1  (4.12) 

then 


G(8)  “  g(8-2n)  when  2n-l<8<2n+l,  n=0,l,2,.. 


.  (4.13) 
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According  to  (3.18),  (3.19),  (3.22),  (3.23),  (4.11)  and  (4.13) 
one  representation  of  the  solution,  which  as  we  w'ill  see  is 
the  one  most  easily  generalized  to  the  nonlinear  case,  is  that 

u  =  g(s)+f(r),  and  c  =  g (s) -f (r) ,  -l<s,r<l,  (4.14) 

where  r(t,X)  and  s(t,X)  are  given  explicitly  by  the  conditions 
that 


t-  (X+2m)  =  r  when  2m-l<t-X<2m+l,  m=0,l,2,..  (4.15) 

while 


t+  (X-l-2n)  =  s  when  2n-l<t+ (X-l)<2n+l ,  n=0,l,2,.. 

(4.16) 

To  determine  f  and  g  in  terms  of  <j>  and  \p  note  that  conditions 
(4.15)  and  (4.16)  imply  that  when  t-0 

m  *  n  =  0,  and  X  =  -r  =  s+1.  (4.17) 

5ince, 

f  *  <p  (X)  and  g  =  i|;(X)  when  t=0,  (4.18) 

f(r)  =  4>(-r)  when  -l<r<0,  (4.19) 

and 

g(s)  =  tj)(l+s)  when  (4.20) 

To  determine  f(r)  over  the  range  0<r<l  and  g(s)  over  the 
range  0<s< 1  we  use  conditions  (4.4)  and  (4.6).  These  yield 

f(r)  =  -g(r-l)  *  -ij/(r)  when  0£r<l,  (4.21) 
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and 


g(s)  =  -f(s-l)  =-<Hl-s)  when  0<s<l.  (4.22) 

In  (4 . 17) - (4 . 22)  we  have  assumed  that  the  boundaries 
X=0  and  X=1  are  not  suddenly  brought  to  rest  at  t=0  so 
that  u  is  single  valued  and  equal  to  zero  at  (X,t)  =  (0,0) 
and  (X,t)  =  (1,0). 

The  statements  (4 . 13) -  (4 . 16)  with  £  and  g  determined 
from  $  and  if/  by  (4 . 19)  -  (4 . 22)  represents  a  complete  solution, 
according  to  linear  theory,  of  the  posed  mixed  initial  and 
boundary  value  problem.  Note  that,  according  to  (4.11)  and 
(4.13),  the  linear  theory  predicts  that  if  the  rates  of 
change  of  F(a)  and  G (8)  are  small,  in  the  sense  (3.2),  over 
the  range  -l£a,8<L  then  they  remain  small  in  this  sense  for 
all  t.  Now,  even  though  the  linear  theory  is  consistent  in 
this  sense,  it  is,  nevertheless,  wrong. ^  As  we  show  in  what 
follows,  even  if  the  rates  of  change  of  F(a)  and  G(8)  over 
the  range  -l<a,8<l  are  small  in  the  sense  (3.2),  ultimately, 
for  sufficiently  large  time,  conditions  (3.2)  are  always 
violated  and  the  finite  rate  theory  must  be  used. 

4.2  Finite  rate  theory. 

Conditions  (3,20)  and  (3.21)  imply  that  in  the  finite  rate 
theory  the  ranges  of  the  characteristic  variables  (a, 8)  are 
(see  fig.  3) 


-IPo-al“  and  "I^o-3-00’ 


(4.23) 


where 


i* 

We  mean  by  this  that  the  linear  theory  of  elasticity 
makes  a  prediction  which  is  in  disagreement  with  the  non¬ 
linear  theory  of  elasticity. 
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P0  =  2  [1-MG (0) ]  and  Q  =  2[1+MF(0)J.  (4.24) 

o 

To  determine  the  signa?  functions  F(a)  and  G(3)  over  this 
range  we  use  an  argument  which  is  similar  to  that  used  for 
the  linear  theory.  The  condition  that  (4.2)  holds  on  X=0 
for  all  t  =  a>0  when,  by  (3.21), 

3  »  a-l+MG (3) ,  =  a-l-MF (a)  (4. 25) 

implies  that 

F (R) +G (R-l -MF (R) )  =  0,  all  R>0.  (4.26) 

Similarly,  the  condition  that  (4.10)  also  holds  on  X=1  for 
all  t  =  3^0  when,  by  (3.20), 

a  «  3-1-MF (a) ,  =  3-l+MG(3)  (4.27) 


implies  that 


F(S-1+MG (S) ) +G (S)  =  0,  all 

To  obtain  a  single  equation  for  F(a) 

a>-ip  we  take 
—  2  o 

R  =  1+MF (R) +S ,  =  1-MG(S)+S,  S > 0  , 
in  (4.26)  to  give 


S>0.  (4.28) 

in  the  range 

R>ip  (4.29) 

— L  O 


F (S+l-MG  (S))+G(S)  =  O',  S>0.  (4.30) 

Conditions  (4.28)  and  (4.30)  state  that  if 

a  =  S-  [1-MG  (S)  ] ,  S>0,  a>~^>0 


(4.31) 
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then 


F (a)+G(S)  =  0  (4.32) 

and 

F (a+2 [1-MG (S) ] ) +G(S)  =  0,  (4.33) 

or,  eliminating  G,  that 

F(a+P(cO)  =  F  (ot)  (4.34) 

when 

P(cc)  =  2 [1+MF (a) ]  ,  all  cc>-b  .  (4.35) 

L  o 

I£,  by  a  procedure  similar  to  that  outlined  above,  F  is 
eliminated,  equations  (4.26)  and  (4.28)  yield  the  equations 

G (8+Q(B) )  =  G($)  (4.36) 


when 


Q(3)  =  2 [ 1-MG (8) ] ,  all  3>-^Q0,  (4.37) 

for  the  variation  in  G. 

Equation  (4.34),  with  P  given  in  terms  of  F  by  (4.35), 
is  a  nonlinear  functional  equation  for  F.  It  states  that 
even  though  F  is  not  periodic  with  fixed  period  2  it  is  in¬ 
variant,  or  repeats,  at  a  sequence  of  times  which  are 
separated  by  an  invariant  time  interval  P=2(1+MF).  In 
particular,  the  zeros  of  F  repeat  at  times  which  are  separated 
by  a  time  interval  2.  When  F(f3)  is  continuous,  the  solution 
to  (4.34)  can  be  expressed  in  terms  of  its  variations  over 
the  range  -yP01^1yPo  w^ere>  according  to  (4.24),  (4.31), 

(4.32)  and  (4.35)" 


P  =  P( 


■— P  ) 
2  oJ 


(4.39) 
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If 


F(cO  =  f(a)  when  "7p0laijP0  (4.40) 

then  the  solution  to  (4.34)  is  given  parametrically  by 

F(a)  «  f  (r) ,  4Po-r-IPo  (4.41) 


where,  when 


(m-i)p0<°t<(m+|)Po>  m=0,l>2 , . .  , 


(4.42) 


r(a)  is  given  implicitly  by 

a  =  r+2m(l+M£(r)) .  (4.43) 

To  check  that  (4. 40) -(4. 43)  does  indeed  represent  the  solution 
of  (4.34)  and  (4.35),  let  am  and  am+1  by  the  values  of  a  in 
the  ranges  (m4)P0<am<(m+^)P0  ctn.*  (m+^)Po<am+1<(m+|)P0  which 
are  the  images  under  the  mapping  defined  by  (4.42)  and  (4.43) 
of  the  same  parameter  r  so  that,  by  (4.40), 

F(<W  =  F(aJ  ■  £W-  (4.44) 

Then,  since 


and 

=  r+2m(l+Mf (r)) 

am+1  =  r+2 (m+1) (1+Mf (r) ) 

> 

y 

(4.45) 

am+l~am  -  2 (1+Mf (r))  =  P(oB)  -  P(«n+1) 

(4.46) 

Equations  (4.44)  and  (4.46)  state  that  F  is  equal  at  two 
values  of  a  whose  difference  in  P.  This  agrees  with  the 
statements  (4.34)  and  (4,35). 

In  contrast  to  the  linear  theory,  the  solution  represented 
by  (4.40) - (4.43)  predicts  that  even  though  the  level  of 
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|MF(a)  1  is  bounded  by  its  level  over  the  range  -^PQ£a<_yP  , 
the  level  o£  |MF'(a)|  is  not  bounded  by  its  level  over  this 
range.  In  fact,  at  some  a=ctc,  which  is  an  image  under  the 
mapping  (4.43)  of  some  point  r=rc  at  which  Mf'(r.J<0, 

|MF'(ctc)|  is  unbounded.  This  implies  that  a  shock  forms  at 
X=0  at  time  t=aclor  somewhere  else  in  the  panel  before  time 
t„.  (Kow  the  analysis  must  be  modified  once  shocks  form  is 
discussed  in  part  II.)  To  see  this  note  that,  by  (4.41) 
and  (4.43), 

F»(o)  =  f'(r)||»  f ' (r) [l+2mMf 1 (r) ] ~1.  (4.47) 

Equation  (4.47)  implies  that  | MF * ( a) |  first  becomes  unbounded, 
compared  with  |Mf*(r)|,  at  some  time 


tc  =  “c  =  tc+2mc(l+Mf (rc)) ,  (4.48) 

where  mc  is  the  least  integer  which  satisfies  the  condition 
that 


-2Mf ' (r  )  =  m-1 

v  QJ  C 


(4.49) 


for  some  rc  in  the  range  -^-P0£rc£jP0 .  Since,  as  we  will  show 
when  we  determine  f(r)  in  terms  of  conditions  in  the  panel 
at  t=0, 


(4.50) 


there  is  always  some  rc  which  satisfies  (4.49)  for  some  mc>0. 

The  signal  function  G(3)  can  be  determined  from  conditions 
(4.36)  and  (4.37)  by  an  argument  similar  to  that  used  to 
determine  F (a) ,  If 

G(g)  -  g (3)  when  ,  (4.51) 


then 


where,  when 


G(3)=g(s),  -|qo<s<|qq 


Cn4)Qo-3-(n+I)Qo’  *=0,1,2,...  ,  (4. S3 


s ( 33  is  given  implicitly  by 


3  =  s+2n(l-Mg(s)) 


(4.54 


If  the  statements  (4.41)-(4.43)  and  (4 . S2) - (4. 54)  are 
used  to  determine  F(o)  and  G(3),  (3 . 18) -  (3 . 21)  can  be  used 
to  represent  the  disturbance  in  the  panel  in  terms  of  f  and 
g.  It  is  given  by 


u  =  g(s)+f  (r) ,  and  c  =  g(s)-f(r),  (4.55 


where 


In 


■fPo-ri?o>  and 


in  the  range 


r(t,X)  is  given  implicitly  by 


(4.56; 


"I5 Pq-* 4PoX-(m+I) Po ’  m=0,l,2,..  ,  (4.57; 


t- (X+2m)  =  r+Mf (r) (X+2m) ; 


(4.58; 


in  the  range 


»  n=0 ,1,2  ..  (4.59] 


s(t,X)  is  given  implicitly  by 


t+(X-l-2n)  =  s+Mg (s) (X-l-2n) .  ■  (4. 601 


3 
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4 . 3  Determination  of  f  and  g  in  terms  of  conditions  in  the 
panel  at  t=0 

If,  when  t=0, 


f  =  <}>(X)  and  g  =  ^(X),  0 <X<  1 , 

then  f(r)  is  given  parametrically  by 

!  cj>(-R)  when  r  =  R[l+M(j)(-R)  ]  ,  -1<R<0, 
f(r)  -  .j 

j  -iKR)  when  r  =  R[l-Mi|>(R)  ]  ,  0 <R< 1 , 


(4.61) 


-  ip  <r<0 ; 

2  o - ’ 

0<r<ip  : 

- 2  o 

(4.62) 


g(s)  is  given  by 

iKl+S)  when  s  =  S [1-MiKl+S) ]  ,  -1<S<0,  -iq  <s<0; 
g(s)  --  2  o-- 


(4.63) 


-<j>(l-S)  when  s  =  S[l+M<f>(l-S)  ]  ,  0<S<1,  0<s4q  . 

v.  “  0 

In  (4 . 61) - (4 . 63)  f  and  g  are  continuous  because  it  is  assumed 
that  the  panel  is  not  suddenly  brought  to  rest  at  t=0  so  that 


iJ>(o)+<j>(o)  =  0,  and  \J>(1)  +  ({>(1)  =  0.  (4.64) 


The  results  of  linear  theory,  described  in  section  (4.1),  are 
recovered  by  formally  taking  M=0  in  (4 . 55) - (4 . 64) . 

The  arguments  used  to  establish  the  results  (4.62)  and 
(4.63)  parallel  those  already  used  in  section  (4.1)  for  the 
linear  theory.  To  obtain  the  first  part  of  (4.62)  first 
note  that,  according  to  (4.57)  and  (4.58),  when  t=0  and  0<X<1 

m  =  0  and  r  =  -X[l+Mf(r)].  (4.65) 

However,  by  (4.61),  when  t=0 

f  (r)  =  <j»(X),  0<X<1,  (4.66) 

which,  together  with  (4.65),  imply  the  first  part  of  (4.62). 
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Similarly,  to  obtain  the  first  part  of  (4.63)  note 
that  (4.59)  and  (4.60)  imply  that  when  t=0  and  0<X<1 , 

n  =  0  and  s  =  (X-l) [1-Mg (s) ] .  (4.67) 

However,  by  (4.61),  when  t=0 

g(s)  =  <KX)  (4.68) 

which,  together  with  (4.67),  imply  the  first  part  of  (4.63). 

It  remains  to  establish  the  second  parts  of  (4.62)  and 

(4.63).  To  determine  f(r)  in  the  r ange  0<r<^-P  we  use  (4.26). 
This  relates  f  at  any  argument  in  this  range  to  g(s)  at  a 
negative  argument  where  it  is  given  by  the  first  part  of 

(4.63).  In  fact,  according  to  (4.26) 

f(r)  =  -g(s)  (4.69) 


when 


r  =  s+l-Mg(s),  0<_r<'jPQ,  -  |qo<s<0.  (4.70) 

Equations  (4.69)  and  (4.70)  together  with  the  representation 

(4.63)  for  g(s)  in  the  range  -^Q^s^O  yields  the  second  part 
of  (4.62).  Finally,  to  determine  g(s)  in  the  range  0_<s<^Qo 
we  use  (4.28)  and  the  first  part  of  (4.62). 

The  statements  (4 . 55) - (4 . 60)  with  f(r)  and  g(s)  given 
by  (4.62)  and  (4.63)  provide  a  complete  algebraic  description 
of  conditions  in  any  small  amplitude  finite  rate  disturbance 
in  terms  of  quite  general  conditions  in  the  panel  at  t=0. 

In  section  6  we  show  that  these  statements  are  still  valid 
away  from  weak  shocks  whose  influence  is  to  attenuate  the 
disturbance . 


4 . 4  Determination  of  f  and  g  for  the  impact  problem. 


A  problem  of  practical  interest  is  the  determination 


-30- 


of  the  motion  produced  in  an  elastic  panel,  or  bar,  which  is 
rigidly  bonded  at  X=1  after  the  other  end,  X=0,  is  suddenly 
moved  and  then  held  fixed.  The  analagous  problem  in  gas 
dynamics  is  the  determination  of  the  flow  in  a  gas  filled 
tube^  which  is  closed  at  one  end}  when  a  piston  is  suddenly 
moved  and  then  brought  to  rest  at  the  other.  Before  the 
front  of  the  disturbance  reaches  the  boundary  X=l,  it  is 
generated  by  the  simple  wave  which  was  described  in  section 
2.  If  the  duration  of  the  motion,  e,  of  the  boundary  X=0 
is  short  in  the  sense  that 


e<2 ,  (4.71) 

that  is,  if  u|0  at  X=0  only  over  a  time  interval  which  is 
less  than  that  taken  by  a  sound  front  to  leave  and  return  to 
X=0  after  reflection  from  X=l,  the  signal  functions  f(r) 
and  g(s)  can  readily  be  expressed  in  terms  of  the  motion 
of  the  boundary  X=0. 

If  at.  X=0 


u  =  U(t) , 


(4.72) 


where 


U(-e)  =  U(o)  =  0 


(4.73) 


and  where 


U=0  outside  the  range  -e<t<0 


(4.74) 


then,  for  -e<t<l-e 


F(o)  =  Ufa)  and  G (3)  5  0  (4.75) 

so,  by  (3 . 18) -  (3 . 21) ,  the  disturbance  is  generated  by  a  simple 
wave.  For  t>_-l  (<l-e)  the  disturbance  is  described  by  (4.55)- 
(4.60)  with 


f(r)  =  i  “  "  i  (4.76) 

|U(a-2)  when  r  =  a+2MU(o-2),  0<a<aQ ,  0<r<|PQ, 

where 


jP0+MU(- jPQ)  =  1  and  a0  =  |po-2;  (4.77) 


with 


g(s)  =  -U(t-1),  when  s  =  t+MU(t-1),  -1<t<1,  -1<s<1. 

(4.78) 


The  statements  (4 . 76) - (4 . 78)  are  readily  established  by  using 
arguments  similar  to  those  used  in  section  (4.3). 
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5.  DISCUSSION  AND  APPLICATIONS 

In  section  3  we  pointed  out  that  the  linear  approximation 
to  F  as  a  function  of  (t,X)  is  valid  only  if  |MF(t) |  and  |MF'(t)| 
are  small  compared  with  unity.  However,  as  for  the  standing 
wave  problem  discussed  in  section  4,  in  most  applications  F 
is  not  known  a  priori  but  must  be  calculated  as  a  solution 
to  some  difference  equation.  An  obvious  question  which  arises 
is  what  are  the  restrictions  which  must  be  imposed  on  F  which 
allow  the  difference  equation  to  be  replaced  by  the  difference 
equation  which  is  predicted  by  linear  theory.  There  is  no 
reason  to  assume  that  they  are  the  same  as  those  which  allow 
F(;.)  to  be  approximated  by  F(t-X):  the  restrictions  will 
depend  on  the  particular  problem  being  investigated.  To 
illustrate  this  point,  we  derive  the  conditions  which  allow 
the  solution  to  the  nonlinear  difference  equation  (4.34)  and 
(4.35),  which  is  given  by  (4 . 41) - (4 .43) ,  to  be  approximated 
by  the  solution  to  the  linear  difference  equation  (4.9). 

For  convenience  we  take  PQ=2,  (this  is  always  possible 
by  a  suitable  choice  of  reference  configuration).  Then,  (4.40) 
and  (4.43)  imply  that 

F(t)  =  f(r),  -l<r<l,  (5.1) 

where  r(t)  is  given  by 

t  =  r+2m(l+Mf (r) ) .  (5.2) 

The  linear  theory  predicts  that 

FL(t)  =  f(rL),  (5.3) 

where 

rT  =  t-2m  =  r+2mMf(r).  (5.4) 

Jj 

According  to  (5.1) -(5. 4),  and  the  mean  value  theorem  of 
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differential  calculus, 


jf(rL)-f(r) 

FltJ- 

2mMf ' (0) | ,  (5.5) 


for  some  0  between  r  and  r  .  According  to  (5.5)  only  when 

L 


max . | 2mMf 1 (6) | <<1 


(5.6) 


does  the  solution  FL(t)  to  the  linear  difference  equation 
approximate  the  solution  F(t)  to  the  nonlinear  difference 
equation.  Since 


MF'(t)  =  Mf ’ (r)/l+2mMf ' (r) ,  (5.7) 


for  all  times  at  which  (5.6)  is  satisfied  the  small  rate 
assumption  | MF  f ( t) | << 1  is  also  satisfied,  so  that  F(a)  can 
be  approximated  by  F(t-X).  However,  the  converse  is  not 
true.  Even  though  |MF'(t)|<<l  over  some  time  interval  the 
restriction  (5.6)  need  not  hold:  even  though  the  distortion 
of  the  signal  carried  by  the  a-wavelets  may  be  neglected  for 
some  time,  (i.e.  F(a)  =  F(t-X)),  nevertheless,  it  does  not 
follow  that  amplitude  dispersion  (nonlinearity)  may  also  be 
neglected  in  calculating  the  signal  carried.  According  to 
(5.6)  and  (5.7)  a  restriction  on  MF'(t)  which  allows  the 
'elution  of  the  nonlinear  difference  equation  to  be  approximated 
by  the  solution  to  the  linear  difference  equation  is  that 


MF' (t) 

(2m) -1-MF ' (t) 


<<1. 


(5.8) 


Chester  [1963]  has  shown  that  small  amplitude  resonant 
vibrations  can  also  be  investigated  by  a  theory  in  which 
amplitude  dispersion  can  be  neglected  in  the  body  of  the 
transmitting  material  but  not  in  the  calculation  of  the 
signal  which  is  carried,  (which  contains  shocks.')  In  a 
future  paper  we  show,  however,  that  the  more  general  finite 
rate  theory  must  be  used  to  investigate  conditions  in  the 
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resonant  band  and  also  to  investigate  the  approach  to  resonance. 

We  have  not  shown  that  the  solution  to  the  nonlinear 
difference  equation  (4.34)  and  (4.35)  does,  in  fact,  itself 
approximate  the  solution  to  the  exact  equation  satisfied  by 
F  in  the  limit  when  MF+0.  Presumably,  this  question  could 
be  resolved  by  using  the  hodograph  representation  of  the  solution. 
However,  in  partial  support  of  our  assumption  that  it  does, 
we  note  that  there  is  an  equation  of  state  for  which  the  difference 
equation  can  be  solved  exactly  so  that  the  question  can  be 
answered.  This  is  the  equation  of  state  for  which  the  co¬ 
efficient  of  F'(ot)  in  (3.9) 

_3_ 

[a(c) ] ~2a 1 (c)  =  constant,  =  M  say.T  (5.9) 

By  a  suitable  choice  of  E  and  M  this  equation  of  state  can 
be  used  to  locally  approximate  many  other  equations  of  state 
when  investigating  small  amplitude  phenomena  which  depend  on 
the  magnitudes  of  both  the  local  tangent  and  curvature  of  the 
stress/strain  curve.  When  (5.9)  is  used,  equation  (3.9) 
integrates  to  give 

(1-^Mc)  {a(a)+MXF  '  (a) }.  (5.10) 

The  function  a(a)  is  determined  as  a  function  of  F(a)  from 
the  conditions  that 

on  X  =  0:  t=a,  u  =  0  and  c  =  -2F(a):  (5.11) 

it  is  given  by 


f 

The  equation  of  state  (S.9)  is  an  excellent  approximation 
to  the  actual  stress/strain  curves  of  a  wide  variety  of  "locking 
materials"  over  the  full  range  of  strain.  It  is  currently 
being  used  by  Cekirge  and  Varley  to  investigate  large  amplitude 
waves  in  such  materials. 
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a  =  [l+MFCa)]"1.  (5.12) 


Similarly, 


||=  (l-JjMc) {S(3)+M(X-1)G'(3)}  (5.13) 

where  "B  is  determined  as  a  function  of  G(3)  from  the  condition 
that 

on  X  =  1:  t  =  3,  u  =  0  and  c  =  2G(3):  (5.14) 

it  is  given  by 


3  =  [1-MG(3)]_1.  (5.15) 

If  now  (5.10)  is  integrated  at  X=1  and  if  the  condition 
that 


on  X  =  1,  c  =  2G (3)  =  -2Ffa)  (5.16) 

is  used  we  obtain  the  condition  that 

on  X  =  1,  t  =  <x+t  +MF(a)[l+^MF(cO]  ,  (5.17) 

for  some  constant  t  .  Similarly,  (5.13)  and  the  conditions 
that 


on  X  =0:  c  =  -2F(o)  =  2G(3)  (5.18) 


implies  that 


on  X  =  0:  t  =  6+t  -MG(B) [l+^MG(0)j  (5.19) 

when  tg  is  a  constant.  If  now  the  three  conditions  (5.11), 
(5.18)  and  (5.19)  on  X=0  are  used,  together  with  the  three 
conditions  (5.14),  (5.16)  and  (5.17)  on  X=l,  then,  by  using 
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arguments  which  are  similar  to  those  used  in  section  4.2,  a 
nonlinear  difference  equation  for  F  is  obtained.  Quite 
remarkably,  if  the  strain  is  measured  so  that  t1+t2=2,  the 
exact  equation  for  F  is  identical  to  that  given  for  F  by 
(4.34)  and  (4.35).  That  is,  equations  (4.34)  and  (4.35), 
which  we  hope  are  good  approximations  for  any  equation  of 
state  in  the  limit  when  MF->0 ,  are  exact  for  the  equation  of 
state  (5.9)  for  any  level  of  MF. 

5 . 1  Associated  pure  initial  value  problem 

One  prediction  of  linear  theory  is  that  the  solution 
to  the  mixed  initial  and  homogeneous  boundary  value  problem, 
which  is  given  in  section  4.1,  can  also  be  derived  as  part 
of  the  solution  to  a  pure  initial  value  problem  in  the 
region  -«><x<«>.  The  initial  values  of  F  and  G  for  this 
associated  problem  are  periodic  functions  of  X  (all  X)  with 
fixed  period  2  and  with 

F  =  f(-X)  for  -1<X<1,  (5.20) 


and  with 


G  =  g(X-l)  for  0<X<2,  (5.21) 

where  f  and  g  are  given  by  (4.62)  and  (4.63)  with  M=0.  The 
form  of  the  dependence  of  f  and  g  on  u  and  c  imply  that  the 
initial  values  of  u  and  c  for  this  associated  problem  are 
also  periodic  in  X  and,  in  addition,  that  u  is  odd,  while  c 
is  even,  with  respect  to  both  points  X=0  and  X=l. 

According  to  the  representation  (4 . 55) -  (4 . 60)  the 
solution  to  the  nonlinear  mixed  initial-boundary  value 
problem  can  also  b>.  constructed  from  the  solution  to  a  non¬ 
linear  pure  initial  value  problem.  Moreover,  the  initial 
values,  4>(X)  and  T(X),  of  F  and  G  for  the  associated  problem 
are  obtained  from  their  prescribed  values,  <f>(X)  and  i|»(X), 

over  the  range  0£X<1  by  the  identical  procedure  used  in 
+■ 

linear  theory.  For,  (4.58)  states  that 
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$(X)=f(r),  -|Po-r-TPo'  (5‘22) 

where  r(X)  is  given  by  the  condition  that 

X+2m  =  -  'f+Mf  ^ry  when  -2m-l<X£l-2m,  m=  ..-1,0,1,.. 

(5.23) 

According  to  (5.22)  4  is  invariant,  or  repeats,  at  values  of 
X  which  are  the  images  under  the  mapping  (5.23)  of  the  same 
parameter  r  but  different  values  of  m.  Any  two  neighboring 
points  of  this  sequence  are  at  a  fixed  distance  of  2  units 
(twice  the  width  of  the  panel)  apart  irrespective  of  the  value 
of  r  so  that 


<D(X+2)  =  $(X).  (5.24) 

In  a  similar  fashion,  it  can  be  shown  that  at  any  later 
time  F  is  still  periodic  in  X  with  fixed  period  2.  It  follows 
that  even  though  the  solution  to  the  associated  initial 
value  problem  is  not  periodic  in  time,  as  in  the  linear 
theory,  it  is  periodic  in  X. 

5 . 2  Distortion  of  a  normal  mode. 

Another  representation  of  any  solution  to  the  linear 
problem  is  as  a  sum  of  its  noninteracting,  nondistorting, 
Fourier  modes.  For  a  nonlinear  system,  however,  the  separate 
Fourier  components  into  which  any  initial  disturbance  can 
be  decomposed  do  interact  so  that  in  the  finite  rate  theory 
this  representation  is  invalid.  In  fact,  even  a  pure  initial 
sinusoidal  disturbance  will  distort  in  time,  generating 
higher  harmonics,  until  shocks  form  .  These  ultimately 
attenuate  the  disturbance.  In  part  II  of  this  report  we 
show  that  the  ultimate  equilibrium  state  of  the  material 
is  a  state  of  rest  in  which  the  strain  is  uniform  and  is 
given  by 
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e  =  f  e(t,X)dX.  (5.25) 

'  o 

(eM  is  independent  of  time  because  the  width  of  a  panel 
which  is  held  fixed  at  both  boundaries  does  not  change.) 

This  final  configuration,  in  which  the  shocks  have  smeared 
out  all  spatial  inhomogeneities,  is  an  appropriate  reference 
configuration  from  which  the  strain  and  stored  elastic  energy 
should  be  measured.  With  this  convention 


(5.26) 


This  reference  configuration  is  attained  as  the  asymptotic 
limit  of  a  sequence  of  deformation  for  which,  between  shocks, 
essentially 


9u  „  9e  . 
9t  “  9X  = 


(5.27) 


while 


9e  _  9u 
9t  9X 


constant . 


(5.28) 


Asymptotically,  for  "large"  time,  at  any  fixed  time,  the 
profile  of  u  against  X  is  a  series  of  N-waves  and  that  of  e 
against  X  a  series  of  square  waves;  at  any  fixed  X,  however, 
the  profile  of  u  against  t  is  a  series  of  square  waves  and 
that  of  e  against  t  a  series  of  N-waves. 

In  figures  4  and  5,  as  an  illustration  of  the  theory, 
we  show  how  the  disturbance  corresponding  to  the  initial 
conditions 

u=0  and  Me  =  ^^(sinTTX—)  (5.29) 

(e  is  measured  to  have  zero  mean)  develops.  We  graph  the 
variations  of  Mu  and  Me  as  functions  of  X  at  times  which 
correspond  to  8  and  14  vibrations  according  to  linear  theory. 
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At  these  times  linear  theory  predicts  u  and  e  are  again 
given  by  (5.29).  Note  how  the  profiles  are  approaching  the 
asymptotic  profiles  given  by  (5.27)  and  (5.18)  and  that 
there  is  a  tendency  to  equipartion  energy  between  kinetic 
and  potential  energy.  This  will  be  discussed  in  greater 
detail  in  part  II. 
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